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By means of detailed Monte Carlo (MC) simulations, we have presented dynamic phase transi¬ 
tion (DPT) properties of ferromagnetic thin-films. Thermal variations of surface, bulk and total 
dynamical order parameters (DOP) for a Him and total order parameter for the films with differ¬ 
ent thicknesses have been examined. Opposite regimes of the critical value of reduced exchange 
interaction (surface to bulk ratio) R c at which the critical temperature becomes independent of 
film thickness L has been also taken into consideration. The average magnetizations of each layer 
is reversed in these regimes. Based on the results, we have confirmed that the system represents 
a crossover behavior in between ordinary to extraordinary transition in the presence of surface 
exchange enhancement. 
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I. INTRODUCTION 

Magnetic properties of free surfaces drastically differ 
from the bulk material, because the free surface breaks 
the translational symmetry (i.e. surface atoms are em¬ 
bedded in an environment of lower symmetry than that 
of the inner atoms and consequently the exchange con¬ 
stants between atoms in the surface region may differ 
from the bulk value). The surface enhancement phe¬ 
nomenon in finite magnetic materials has attracted con¬ 
siderable amount of interest for both experimentalists [l|- 
[f| and theorists EM- 

Applied oscillating magnetic fields, depending on the 
competition between the two time scales, namely the os¬ 
cillation period P of the external perturbation and the 
relaxation time r of the sample, a dynamic symmetry 
breaking may take place causing a DPT. There are two 
cases due to the competition between these time scales: 
P < t and P > t. In the first case, the system can¬ 
not relax within a complete cycle of the magnetic field 
oscillation, hence the instantaneous magnetization M{t) 
oscillates in time around a nonzero value corresponding 
to dynamically ordered state (i.e. dynamic ferromagnetic 
phase). In other case, M(t) can follow the external field 
with some delay, and the system exhibits a dynamic para¬ 
magnetic behavior. The relaxation time r can be con¬ 
trolled by supplied energy with several different ways: 
The agency of an adjustable parameter such as the field 
amplitude, strength itself, the type of the exchange in¬ 
teractions, and the temperature. The DPT point can 
be controlled by tuning mentioned competing factors to¬ 
gether with the time period of external field. 

Experimental point of view, Schierle and coworkers ob¬ 
served that the magnetizations of the outermost layers 
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in EuTe(lll) films decrease significantly differently from 
those of bulk layers [I| . Violbarbosa and coworkers found 
that the formation of the blocks of layers with robust 
magnetic structure whereas the interblock interactions 
are relatively weak in fcc-Fe on Cu(001) film Q. More¬ 
over, the enhanced surface magnetism has been the focus 
of such systems. For instance, Gd film has been inves¬ 
tigated experimentally. The thickness-dependent spin- 
polarized electronic structure of strained ultrathin and 
thin films of Gd has been investigated by Waldfried et 
al. E3- They found that the surface magnetic struc¬ 
ture dominates the magnetic ordering of the ultrathin Gd 
films. With decreasing thickness some bulk bands exhibit 
increasingly more passive magnetic behavior. Skomski 
and coworkers also found that Gd films exhibit a mag¬ 
netic surface transition which occurs at about above the 
bulk Curie temperature [H| . 

In sense of dynamic phase transitions, a great deal of 
theoretical efforts has also been devoted to the investiga¬ 
tion on such systems. The details of surface enhancement 
phenomenon for the films were subjected to an exter¬ 
nal oscillatory field have been intensively propounde d by 
Akta§ et al. by using effective field theory (EFT) flff . 
The general trend of frequency dispersion belongs to crit¬ 
ical temperature coordinate of the special point for dif¬ 
ferent frequency and amplitude values has been demon¬ 
strated in their work. Nonequilibrium phase transition in 
the kinetic Ising model on a two-layer square lattice has 
been examined by Canko et al. f2p| . Dynamic phase dia¬ 
grams have been constructed in the plane of the reduced 
temperature versus the amplitude. Similarly, dynamic 
magnetic behavior of a mixed Ising system on a bilayer 
square lattice has been investigated by Erta§ and Keskin 
121| . They presented the dynamic phase diagrams in the 
reduced temperature and magnetic field amplitude plane 
and the effects of interlayer coupling interaction on the 
critical behavior of the system have been investigated in 
their work. 
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In recent series of works by Pleimling and coworkers, 
surface criticality at a DPT and surface phase diagram 
of the three-dimensional kinetic Ising model has been 
elucidated. In the first one of these studies, Park and 
Pleimling found that the nonequilibrium surface expo¬ 
nents do not coincide with those of the equilibrium crit¬ 
ical surface (22| . In addition, in three space dimensions, 
the surface phase diagram of the nonequilibrium system 
differs markedly from that of the equilibrium system. 
The values of the critical exponents have been deter¬ 
mined through finite-size scalin g by Park and Pleimling 
in their followup investigation J23|. Their results have 
showed that the studied nonequilibrium phase transition 
belongs to the universality class of the equilibrium three- 
dimensional Ising model. The surface phase diagram of 
the three-dimensional kinetic Ising model below the equi¬ 
librium critical point subjected to a periodically oscillat¬ 
ing magnetic field has also been presented by Taucher 
and Pleimling {24[- They presented that surface phase 
diagram that in parts strongly resembles the correspond¬ 
ing equilibrium phase diagram, with an ordinary transi¬ 
tion, an extraordinary transition, and a surface transi¬ 
tion. These three lines meet at a special transition point. 
For weak surface couplings, however, the surface does not 
order. 

In this regard, our task in the present work is to shed 
some light on the DPT properties, -especially the evolu¬ 
tion of crossover point with field amplitude- of ferromag¬ 
netic thin-films in the presence of ac driving fields. In 
this present paper, the DPT properties in the presence 
of external oscillatory field of the system are studied by 
MC simulation. The layout of the work is as follows: 
Section 2 describes the model and the MC simulation 
scheme, the numerical results are reported in Section 3, 
the paper ends with concluding remarks in Section 4. 


II. SIMULATION 

We consider a ferromagnetic thin film with thickness 
L described by spin-1/2 Hamiltonian 

U = ~ h ^) Si C 1 ) 

<«> * 


In order to simulate the system, we employ the 
Metropolis MC simulation algorithm [25], [13 to Eq. ([1]) 
on an N x TV x L simple cubic lattice where N = 70 and we 
apply periodic (free) boundary conditions in direction(s) 
parallel (perpendicular) to film plane. We have studied 
ultrathin-films with thickness L = 3,4,5 together with 
a relatively thicker thin-film L = 20 to observe average 
magnetizations of each layer for selected some system 
parameters. For simplicity, the exchange couplings are 
restricted to the ferromagnetic case. 

Configurations were generated by selecting the sites 
in sequence through the lattice and making single-spin- 
flip attempts, which were accepted or rejected according 
to the Metropolis algorithm, and N x N x L sites are 
visited at each time step (a time step is defined as an 
MC step per site or simply MCS). Data were generated 
over 50 independent sample realizations by running the 
simulations for 50000 MCS per site after discarding the 
first 25000 steps. This amount of transient steps is found 
to be sufficient for thermalization for the whole range of 
the parameter sets. Throughout the analysis, oscillation 
period of the external field is kept fixed as P = 100. 

Our program calculates the instantaneous values of the 
bulk and surface magnetizations M s and M b , and the 
total magnetization Mr at time t. These quantities are 
defined as 


1 jv„ JV & 

M s{t) = = TT 5Z S P 


N s 
M T (t ) = 


N b . 

i—i j=i 

N s M s (t) + N b M b (t) 

N s + N b 


( 2 ) 


where N s and N b denote the number of spins in the sur¬ 
face and bulk layers, respectively. From the instanta¬ 
neous magnetizations, we obtain the related order pa¬ 
rameters as follows |27l]: 


Q s 


1 

P 


M s (t)dt , Q b = p 
Qt = -p f> M T {t)dt 


M b (t)dt, 


( 3 ) 


Using Eq. ©, we calculate the total energy per spin 


where s,; = ±1 is a two-state spin variable, and Jij is the 
nearest neighbor interaction energy. The summation in 
the first term is taken over only the nearest neighbor in¬ 
teractions whereas the summation in the second term is 
carried out over the all lattice sites. In the second term, 
h(t) = ho sin(o;t) represents the oscillating magnetic field, 
where ho and w are the amplitude and the angular fre¬ 
quency of the applied field, respectively. The period of 
the oscillating magnetic field is given by P = If 

the lattice sites i and j belong to one of the two sur¬ 
faces of the film we have Jij = J s , otherwise Jij = J b , 
where J s and J b denote the ferromagnetic surface and 
bulk exchange interactions, respectively. 


Etot ~ P{N s + N b )j Udt (4) 

Consequently, the specific heat is defined as 

C = ^F' (5) 

We also note that the value of the bulk exchange interac¬ 
tion J b is fixed to unity, and we also use the normalized 
surface to bulk ratio of exchange interactions R = J s / J bl 
as well as the reduced field amplitude Hq = ho/J b , and 
reduced temperature 0 = ksT/ J b . 
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III. RESULTS AND DISCUSSION 

Based on the results, we mainly focus on the effect of 
external oscillatory field amplitude in surface enhance¬ 
ment phenomenon. First, in order to obtain a general 
insight on DPT characteristics, we plot Fig. © and © 
respectively. In Fig. © the thermal variations of sur¬ 
face, bulk and total order parameters for the films with 
thickness L = 3 are shown. Following this, total dy¬ 
namical order parameters for the films with three differ¬ 
ent thicknesses are shown in Fig. ©■ We restrict our 
discussions for two value of reduced exchange R = 0.25 
and 2.75. Both in Fig. © and ©, the transition point 
increases with reduced exchange interaction at constant 
values of the other system parameters. Hence, relatively 
more thermal agitation is needed to make the system dy¬ 
namically disordered for more interaction. Moreover, at 
a constant temperature, surfaces are weakly ordered due 
to the scarcity of dipole-dipole interaction per site for 
R = 0.25 value. For R = 2.75, surfaces dominate against 
bulk due to the surface enhancement. From Fig. ©, 
we see that both surfaces and bulk layers of magnetic 
thin-films exhibit a phase transition at a certain critical 
temperature independently from the value of R. 



neighbor bonds, and a DPT can be observed at low tem¬ 
peratures. When we fixed rest of the parameters except 
reduced exchange and compared the strength of the or¬ 
der parameters for the films with different thicknesses in 
two different value of R (namely R = 0.25 and 2.75), we 
see the hierarchically sequence reversal in Fig. © (from 
(a) to (d) and (b) to (e) and (c) to (f)). In between these 
two value, there should be a critical point of the reduced 
exchange at which all the layers seem to oscillate in phase 
independently from the thickness. 




0 
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FIG. 1. (Color online) Surface, bulk and total order parame¬ 
ters for a film L = 3 layers for two different regime of reduced 
exchange interaction. 

From Fig. ©, one can easily see the effect of external 
field amplitude Hq for fixed value of the system parame¬ 
ters from the panel (a) to (c) and (d) to (f). Critical tem¬ 
perature exhibits a decreasing behavior with increasing 
Hq as a consequence of the well-known following physical 
mechanism: For small amplitude values, the energy sup¬ 
plied by the external oscillatory field cannot break the 
ferromagnetic energy induced order due to the nearest- 
neighbor exchange coupling through the system at low 
temperatures. Hence, a DPT cannot be observed un¬ 
less a relatively large amount of thermal energy is sup¬ 
plied to the system. As the field amplitude increases, 
it becomes dominant against the ferromagnetic nearest- 


FIG. 2. (Color online) Average magnetization in static case 
(a) and (d), and the total order parameters for the films with 
three different thicknesses L = 3, 4, 5. and for three different 
values of Ho- 

In order to make the aforementioned phenomenon 
more clear, we plot the dynamical order parameters of 
each layers for a film with L = 20 in Fig. ©. For 
this purpose, we choose a constant temperature value at 
which the system is well-below the transition point and 
the thermal fluctuations can be ruled out. The magneti¬ 
zation M(t) cannot follow the external field h(t) (r > P 
case) for each selected field amplitude values H 0 , conse¬ 
quently the dynamic ferromagnetism is enhanced. Re¬ 
duced exchange varies from R = 1.0 to 2.0 including the 
critical value of itself. So, we have qualitatively different 
two regimes: R < R c and R > R c . Below R c the in- 
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FIG. 3. (Color online) Average magnetizations profiles for 
a film with L = 20 layers. The number accompanying each 
curve denotes several values of reduced exchange varies from 
1.0 to 2.0 including R c . Each layer labelled by k. 

ner layers are highly ordered compare to surfaces. This 
can be briefly explained as follows: In middle of the film, 
there are relatively more neighboring per magnetic sites 
which causes locally larger magnetic interaction. So it 
becomes more difficult to follow the external field for any 
spin. Surface spins are embedded in an environment of 
lower symmetry than that of the inner atoms. The ex¬ 
change constant between atoms in the surface region may 
differ from the bulk one. This regime corresponds to a 
surface type of magnetic ordering. The opposite of the 
above scenarios can be considered also. Above R c , in 
the inner layers, although there are more neighboring, 
there are far fewer exchange constant per magnetic sites 
which causes relatively smaller magnetic interaction than 
that of the surface one. R c plays the main role to obtain 
the frontier of this crossover. The free surface cannot 
break the translational symmetry since magnetic prop¬ 
erties of the free surfaces exactly overlap with the bulk 
one at R c . We can say more generally that the deficiency 
of the interaction per surface spin can be compensated 
by increasing the modified exchange interaction strength. 
Moreover, the effect of external field can be also seen by 
following the panels from (a) to (c). 

In obtaining the critical frontiers depicted in 


( ksT c / J — R) plane, we evaluated the thermal variation 
of specific heat for a given set of system parameters. A 
typical example is shown in Fig. (jTj for the films with dif¬ 
ferent thicknesses as L = 3,4, 5. The temperature values 
corresponding to the maxima of specific heat curves are 
the transition temperatures. From the panel (a) to (c) 
and (d) to (f), field amplitude H 0 changes. Also, from 
(a) to (d), (b) to (e) and (c) to (f) reduced exchange 
has two different value. DPT points shift towards the 
lower temperature with increasing field values as well 
as they shift towards higher temperature with increas¬ 
ing reduced exchange. The related detailed story has 
been explained above. Similarly, the effect of reduced 
exchange on specific heat peaks is easy to understand: 
The stronger dipole-dipole interaction make more contri¬ 
bution to total energy. Hence, more thermal agitation is 
needed to make the system dynamically disordered. In 
Fig. ©, the crossover behavior can be seen easily when 
R changed from R < R c to R > R c (namely, from (a) 
to (d), (b) to (e) and (c) to (f)). Below R Cl thicker film 
has more nearest-neighbor interaction per site, this cre¬ 
ates more contribution to energy. Consequently, both 
strength of the peak and corresponding critical temper¬ 
ature are relatively higher than the others. 

In order to obtain a general overview of the non¬ 
equilibrium phase diagram in (ksT c / Jb — R) planes. For 
this purpose, in Fig. ([5]) we plot the critical temperature 
versus R with selected film thickness values L = 3,4,5 
and for three selected values of field amplitude H 0 . Since, 
the temperature values at which the specific heat curves 
exhibit a sharp maximum correspond to the transition 
temperature of thin film, critical temperature values have 
been obtained by examining the thermal variation of spe¬ 
cific heat curves (a selected set has been given in Fig. 
©)■ Fig. © represents a characteristic phenomenon 
peculiar to thin film systems. Namely, due to the exis¬ 
tence of reduced surfaces, there exists a special value of 
surface to bulk ratio of exchange interactions Rc at which 
the transition temperature of the film becomes indepen¬ 
dent of thickness L. Simply we can say that the curves 
with different film thicknesses intersect each other. This 
fully supports a recent study fl9l ] where in the framework 
of an EFT the existence of a special transition point was 
predicted. The critical temperature value of crossover 
point for Ho = 0.0 (static case) is in a good agreement 
with previous studies [HI, HD, HH • However, variation of 
R c as a function of H 0 is very slow according to Fig. m, 
and we see that the location of R c barely deviates from 
its equilibrium value with increasing Hq. This deviation 
has been reported before by Yiiksel [l2|- The discussion 
on existence of this kind of deviation is an academic is¬ 
sue and this may be due to insufficient data and cannot 
be located accurately with such an approach with less 
effort. This was also reported in a MC simulation treat¬ 
ment of surface critical phenomena by Hasenbusch [20l |. 
From panel (a) to (c) the effect of field amplitude Hq on 
the transition characteristics of the film can be seen and 
it is also straightforward as stated before: Greater the 
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FIG. 4. (Color online) The specific heat versus temperature 
for the films with different L = 3,4, 5 thicknesses. 


FIG. 5. (Color online) The critical temperature versus re¬ 
duced exchange interaction in phase planes for the films with 
different L = 3,4, 5 thicknesses. 


amplitude Hq means more energy transferred to the sys¬ 
tem over half cycle by the external oscillatory field and 
this makes transition from dynamically ordered to disor¬ 
dered phase more easy. For R < R c , we have ordinary 
transition behavior where the bulk magnetism is domi¬ 
nant against the surface magnetism whereas for R > R Cl 
the surface may exhibit enhanced magnetic behavior in 
comparison with bulk. This is called extraordinary tran¬ 
sition. Moreover, as shown in Fig. ([5j), for R < R Cl 
thicker films have greater transition temperatures while 
for R > R c , the transition temperature of the film de¬ 
creases with increasing thickness. The results indicate 
that the well-known surface enhancement properties of 
the system may change its characteristics in the presence 
of an external oscillatory field. 


IV. CONCLUSION 

In conclusion, we have applied MC simulations to 
study the DPT characteristics in thin ferromagnetic films 
in the presence of oscillating magnetic fields. The fore¬ 
most results obtained from simulation data can be sum¬ 


marized as follows: We first investigate the thermal vari¬ 
ations of related order parameters for the films with dif¬ 
ferent thicknesses. The effect of the field amplitude and 
reduced exchange on the previous standard arguments 
for static case has been propounded. In the vicinity 
of the dynamic ferromagnetic-paramagnetic phase tran¬ 
sition temperature, specific heat curves exhibit a sharp 
peak which becomes more apparent for sufficiently high 
reduced exchange values (R > R c regime) values. The 
thinner films in the absence of enhanced surfaces (R < R c 
regime) with high field amplitudes exhibit a weak peak 
in relatively lower temperatures. 

According to our findings, an increment of the field 
amplitude causes a decreasing in corresponding temper¬ 
ature coordinate of the crossover point in (ksT c / Jb — R) 
planes. Critical value of surface to bulk ratio of exchange 
interactions R c at which the transition temperature is in¬ 
dependent of film thickness is not apparently responsive 
to varying field amplitude values, but exhibits slow vari¬ 
ation as a function of Hq. We confirmed the general 
trend was generated by using EFT calculations before 
0 - Hence, we can say that the evolution of a crossover 
is not from the limitation of EFT. 

We hope that this study will shed light on further in¬ 
vestigations of the dynamic nature of critical phenomena 
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in pure crystalline ferromagnetic thin films and will be 
beneficial from both theoretical and experimental points 
of view. 
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